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A simple method to obtain a canonical partition function 
for one dimensional lattice gas model is presented. The sim- 
plification is based upon rewriting a sum over all possible 
configurations to a sum over all possible numbers of clusters 
in the system. 



I. CANONICAL ENSEMBLE 

The behavior of canonical ensemble is fully described 
by its canonical partition function. In the case of the 
lattice gas model the partition function can be written in 
the form 



Z(L,N,/3) 



J2 exp(-£[{ e! :}]/3), (1) 



N 



where L is the number of sites in the system, TV is the 
number of particles, (3 = (UbT)^ 1 , ks being the Boltz- 
man constant, -E[{ei}] is potential energy of the system, 
e, = 1 if the i-th site is occupied and = otherwise. 
The main problem when calculating the canonical par- 
tition function is accounting only for states with total 
number of particles N, that is to calculate degeneracies 
of the energy levels. 

Here I present a method for obtaining the parti- 
tion function as well as clustering probability for one- 
dimensional (ID) lattice gas model, assuming that the 
interaction between molecules is only of nearest-neighbor 
type and that the binding energy s s is constant for each 
site except for the leftmost site i = 1 and the rightmost 
site i — L where the binding energies are equal to e s + El 
and e s + sr, respectively. 

First, it is useful to realize, that in the case of nearest 
neighbor interaction, the interaction energy depends only 
on the number of clusters, and not on their positions or 
their lengths. Then, the problem of finding Z(L, N, (3) 
can be formulated as a problem of finding all possible 
configurations of k clusters. 

This is a unique property of one-dimensional systems 
and is based on the fact that the boundary of a ID 
cluster always consists of two points. Two- and three- 
dimensional clusters do not posses this property; the 
boundary of two or three-dimensional clusters depend not 
only on the cluster size but also on the configuration of 
particles within the cluster, i.e., the cluster shape. 

In the following it will be assumed that k > 1, L > 
N > 1. It is suitable for further calculations to intro- 
duce function F{k,n), which counts the number of ways 



in which k integer numbers a.; > 1, i — 1, . . . , k can be 
chosen such that Ei=i a * = n - Thus, it is possible to 
write 



i-k+l 



F(k,n) 



J2 1= F{k-l,n-m). (2) 



Ek 
i=l a i = n 

This reccurent formula gives 
F(k,n) = 




if k > 1 A n > k 
otherwise. 



(3) 



Thus, the number of ways to arrange N particles into 
k clusters is given by F (k, N) and the number of ways 
to arrange (L — N) empty sites into kl clusters is given 
by F(k! ', L — N). So that the number of configurations of 
the system with N occupied sites arranged into k clusters 
and (L — N) epty sites arranged into k' clusters can be 
expressed as F(k, N)F(k', L ~ N). 

Now it is suitable to express the canonical partition 
function Z(L, N, (3) as a sum over number of clusters k = 
!,-••) k max , k max being equal to min(L - N + 1, N) 



Z(L,N,f3)= J2z cl (L,N,k,(3)eM(N-k)sP), (4) 

k=l 

where Z c i(L, N,k, (3) is the partition function for a sys- 
tem of length L with N particles in k clusters and k max 
gives the maximum number of clusters the system can 
contain. 



A. System with interaction between particles and 
boundaries 



To calculate the canonical partition function of the sys- 
tem Z(L, N, (3) in the case of interaction of particles with 
the boundaries of the sytem, it is convenient to write the 
function Z c i(L, N, k, (3) as a sum of three parts 

2 

Z cl (L, N, k,0) = ^2 Wi(L, N, k) exp{S t (3), (5) 



i=0 



where 



£2 — £r + £l, 

£1 = /3 _1 log(exp(efl/3) + exp(e L f3)), 

£0 = 0, (6) 
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and £r{£l) is nearest neighbor interaction between right 
(left) end of the system and particle, and Wi{L, N,k),i = 
0, 1, 2 are the numbers of realizations of system with the 
length L containing N particles arranged into k clusters 
and with i clusters attached to the ends of the system. 

In the case when sites 1 and L are occupied (two clus- 
ters attached to the ends of the system, i = 2), we have 
k' = (ft — 1) clusters of (L — N) empty sites and ft clusters 
of N particles so that we can write 

W 2 (L, N, k) = 5l,n + F(k, N)F(k - 1, L - N), (7) 

where 



1 if n = m 
otherwise. 



(8) 



When only one of the sites 1, L is occupied (i = 1), then 
(L — N) empty sites are arranged into k' — k clusters and 
again we have k clusters of N particles and we can write 



Wi(L, N, k) = F{k, N)F(k, L - N). 



(9) 



When neither one of the sites 1,L (i — 0) is occupied 
we are dealing with k! — (ft+1) clusters of (L — N) empty 
sites and k clusters of N particles, so that it is possible 
to write 



W (L, N, k) = F(k, N)F(k + 1, L - N). 



(10) 



Thus, using the equation (Q) together with the equa- 
tions (H|t]|]) and ( [To| ) we can calculate the partition func- 
tion of the system together with all thermodynamics 
quantities like mean energy of the system E(L, N,T), 
thermal capacity C(L, N,T), pressure p(L,N,T) and 
probability P c i(L, N, (3, ft) to find k clusters in the sys- 
tem, which can be written in the form 



P cl (L,N,P,k) 



Zci(L,N,k,0) 
Z(L,N,f3) 



exp ((N - k)ep) . (11) 



To calculate the probability P p (N, L, (3,x) to find 
a cluster of length x it is needed to start with 
W P (N, L, ft, ft', a:), which is the number of clusters of the 
length x when taking into account all possible configura- 
tions of system of the length L with TV particles arranged 
into ft clusters and with (L—N) empty sites arranged into 
ft' clusters. It can be showed that for ft > 2 

W P (N, L, ft, ft', x) = kF(k - 1, N - x)F(k', L - N), (12) 

and for ft = 1 we get 

W p (N,L,l,k',x)=5 N , x (F(k',L-N) + 6 L , N ). (13) 

Using the same arguments and method that were 
used to calculate Z(L,N,f3), it is possible to express 
P p (L,N,j3,x) as 



P p (L,N,P,x) 



1 



Z p (L,N,f3) 
Z c i(L,N,l,0)6 N , x exp(-s/3) 



F(k — 1 N — t) 
+ ^ kZ cl (L,N,k,0) [ ' AA ^exp(-fce/?) 

k=2 



,N) 



(14) 



It is important to point out that Z P (L, N, (3) 
and Z(L,N,(3) are different functions: Z p (L,N,f3) 
is normalization factor needed to satisfy condition 

T,x=i P p( L i N iPi x ) = 1 ' while Z ( L , N ,P) is the Pe- 
tition function for the system of L sites containing N 

particles. 



B. Periodic boundary conditions 

The calculation of the canonical partition function 
Z p bc(L, N, (3) in the case of the periodic boundary condi- 
tions is now very straightforward, it is enough to realize 
that this problem is equivalent to the case of ends in- 
teracting with particles with interaction energies er — 
£l = e/2 when both first and last site are occupied and 
Er = El = if one of these sites is empty. It immediately 
yields result 



Z pbc (L,N,(3) = 53 Z^ bc) (L,N,k,(3)exp((N-k)E(3), 

(15) 



fc=i 



where 



Zf c> (L, N, ft, (3) = S L , N exp(e/3) + ^F(k, L - N)F(k, N). 

(16) 

Due to the periodic boundary conditions the maximum 
number of clusters k ma x the system can contain is now 
equal to min(L — N, N). 

The expression [l6] can be easily understood. The first 
part Ji t jv exp(e/3) is trivial and it represents a special 
case when L = N, i. e. every site occupied. If L > N, 
then due to the periodic boundary conditions the number 
of clusters of particles is equal to the number of clusters 
of empty sites. Thus, ^(ft, N) counts the number of ways 
to arrange TV particles into ft clusters and F(k, L — N) 
does the same with L — N empty sites. Factor L arises 
from the fact, that there are L ways how to position each 
arrangement of clusters, and factor 1/ft is there because 
in the system with the periodic boundary conditions we 
do not know which cluster is first - and there are exactly 
ft ways to choose the first. 

Calculating the probability P^ bc '(L, N, /3, ft) to find ft 
clusters in the system is now easy 
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P$ bc \L,N,[3,k) 



Z pbc {L,N, 0) 



exp((iV-fc)e/3), 
(17) 



and the probability Pp pbc \L , N ', j3 ', x) to find cluster of 
length a: can be written in the form 

Z^ hc) (L,7V,/3) 

(18) 

where Z p phc \L,N,(3) is normalization factor such that 

J2x=i p p PbC \ L > N ,P> x ) = 1 and Wi.j = 0,1,2 can be 
expressed as 

W (pbc) (L, iV, /3, a?) = LQ{L -N- 1)S X , N + exp(E)6 x , L , 

(19) 



wf pbc) (L,N,p,x)= Y^kF(k-l,N-x) 

x [F(k + 1,L-N) + F{k, L - N)] exp(-fe/3), (20) 
and 

W { 2 pbc) {L, N, (3, x) — ^ F(k — 1,L — N) 

fe=3 

x [(k - 2)F(k -l,N-x) + (x- l)F(k - 2, N - x)] 
xexp(-fce/3). (21) 



II. GRANDCANONICAL ENSEMBLE 

It is now possible to calculate properties of the 
open system (grand-canonical ensemble) since the grand- 
canonical partition function can be writen as a sum of 
canonical partition functions over number of particles 

L 

Q(L, 0, M ) = 1 + £ Z(L, N, (3) exp(N(e s + fifi) 

N=l 

where e s is an interaction energy between particle and 
site (in our ID system) different from 1,L. However, in 
the case of the periodic boundary conditions or in the 
case of the noninteracting boundaries it is simpler to use 
the transfer matrix method [jjj . 



particles was presented. Numerical values can be easily 
obtained, since the summations over all possible config- 
urations of N particles on L sites were rewritten to the 
summations over k = 1, . . . , min(i — N + 1, N), where 
min(i — N + 1, N) is the maximum number of clusters in 
the system, this number being considerably smaller than 
the number of all possible configurations of the system. 

The method employs the fact that the boundary of 
each cluster consists of two points. This is a unique 
property of one dimension, so that it is not possible to 
generalize this method to higher dimensions. 

The results can be used in the study of equilibrium 
properties of quasi-lD clusters. An ideal system on which 
it is possible to study low-dimensional physics is provided 
by single- walled carbon tubes (SWNT's). These can be 
filled with various molecules and atoms, and if the radius 
of a SWNT's is sufficiently small, the motion of trapped 
particles is confined along the tube axis. What makes 
SWNT's special is the possibility to open and close the 
ends of nanotubes 0] , providing the opportunity to study 
both open and closed systems. Recently, the presented 
method was succesfully used by Hodak and Girifalco 
for a theoretical study of clustering Cgo molecules encap- 
sulated in carbon nanotubes. This type of system, also 
called peapod, was experimentaly observed by Smith et 
al. |. 

Other interesting systems, which can be used to probe 
the physics of low-dimensional structures and for which 
the obtained results may be applicable are provided by 
the interstitial channels in ropes of SWNT's f| |6| a nd by 
the grooves on the external surface of the ropes 0] , both 
of which can provide high-energy adsorption sites. 

The presented method can also be useful in analysis 
of the stability of one-dimensional metals, which can be 
obtained for example by alkali-metal absorption on semi- 
conductor surfaces || , and in the study of the behaviour 
of classical gases confined in nanopore materials and 
zeolite materials H. 
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III. CONCLUSIONS 

A simple method for calculation of the properties of 
open and closed ID systems with the periodic bound- 
ary conditions and with the boundaries interacting with 



[1] Rodney J. Baxter, Exactly solved models in statistical me- 
chanics, New York, Academic Press, 1982. 

[2] A. Kuznetsova et al, Chem. Phys. Lett. 321, 292 (2000). 

[3] M. Hodak and L. A. Girifalco, to be published in Phys. 
Rev. B. 



3 



[4] B. W. Smith, M. Monthioux, and D. E. Luzzi, Nature 396, 
3239 (1998). 

[5] G. Stan and M. W. Cole, Surf. Sci. 395, 280 (1998). 
[6] W. Teizer et al, Phys. Rev. Lett. 82, 5305 (1999). 
[7] K. A. Williams and P. C. Eklund, Chem. Phys. Lett. 320, 
352 (2000). 

[8] S. Modesti et al, Surf. Sci. 447, 133 (2000). 
[9] T. W. Ebbesen, Phus. Today 49, 26 (1996). 



4 



